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Abstract 



(N 
> 

, The leading and next-to-leading radiative corrections to deep inelastic events with 

I tagged photons are calculated analytically. Comparisons with previous results and 

■ numerical estimations are presented for the experimental conditions at HERA. 

00 

^ ■ 1 Introduction 

Q^l It is well known that the radiative corrections to deep inelastic electron proton scattering due 

to hard real photon emission are very important in certain regions of the HERA kinematic 
domain. In fact, the initial-state collinear radiation leads to a reduction of the projectile 
electron energy and therefore to a shift of the effective Bjorken variables in the hard scattering 
^ i process as compared to those determined from the actual measurement of the scattered 

' electron alone. Therefore, radiative events 

eiPi) + piP) ^ eip2) + i{k) + X + (7) (1) 

are to be carefully taken into account [|T], ^ 0. 

On the other hand, measuring the energy of the photons emitted very close to the incident 
electron beam direction 0, ^, |^, |^ permits to overlap the kinematical region of photopro- 
duction (Q^ k, 0) and the DIS region with small transferred momenta (about a few GeV^) 
within the high energy HERA experiments. Furthermore, these radiative events may be used 
to independently determine the proton structure functions F2 and Fi (and therefore F^) in 
a single run without lowering the beam energies [^, |^ . Preliminary results of an F2 analysis 
using such radiative events were recently presented by the HI collaboration ||^. 

Our aim is to calculate the radiative corrections to neutral current deep inelastic events 
with simultaneous (exclusive) detection of a hard photon emitted very close to the direction 

*Supported by Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologie (BMBF), Ger- 
many. 
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of the incoming electron beam {9^ = p^k < 9q ^ ■ 10"*^ rad). In the case of the HERA 
colhder, the experimental detection of photons emitted in this very forward direction is 
actually possible due to the presence of photon detectors (PD) that are part of the luminosity 
monitoring system of ZEUS and HI. 

Let us briefly review the kinematics for the process under consideration. As the opening 
angle of the forward photon detector is very small, and since we will only consider cross 
sections where the tagged photon is integrated over the solid angle covered by this photon 
detector, we can parameterize these radiative events using the standard Bjorken variables x 
and that are determined from the measurement of the scattered electron, 

2P ■ (pi - P2) 



q2 = 2pi ■ p2 = xyV, with V = 2P ■ pi , (2) 



and the energy fraction z of the electron after initial state radiation of a collinear photon 

2P -{pi-k) e - k° 



V 



(3) 



where e is the initial electron energy, and is the energy seen in the forward photon detector. 

An alternative set of kinematic variables that is especially adapted to the case of collinear 
radiation, is given by the shifted Bjorken variables [^] 

Q =-[Pi-P2-k) , x=—— — , y = — — — . 4 

2P ■{pi-p2-k) P ■ {pi - k) 

The relations between the shifted and the standard Bjorken variables read 

Q = zQ , x = — , y = . (5) 

z + y — 1 z 

The cross-section under consideration in the Born approximation, integrated over the 
solid angle of the photon detector (0 < 6*^ < ^^o, ^0 ^ 1) then takes the following form: 

Z d^CTBorn _ 1 d^(TBorn _ _^p/^ ^ ^g^ 

y dx dy dz y dx dy dz 2it ' 

where 

27ra2(_Q2) r .o.oM^ ( .oM'^\ f 



Q'^xy'^ 



Q2 V Q^Jl + R 



P Lo = Lo - , R = R{x, g2 = 1 + ^ , -1 

l-z 1-z V J 2xFi{x,Q^) 



a{-Q-') = ^-^^^y Lo = ln(^-^j, = 2zpr ■ p, = 2ze'Y {I - c) 
£2 , , Ep{l+(3p) 



Y = — = l-y + xy — , c = cos(piP2) 



2 



^ zeY{l-c) ^ k-M-lE- 

. ^ 2P-{zp^-p2) ^ z{l+Pp)-Y{l+Ppc) 

^ zV z{l + Pp) • ^ ^ 

The quantities F2 and Fi are the proton structure functions; M and m are the proton 
and electron masses, respectively. In the cross-section (|^) we take into account terms pro- 
portional to M^/Q^ , which may be important at low Q^. Note that the neglect of Z-boson 
exchange and 7-Z interference is a good approximation, because we are interested mostly 
in events with small momentum transfer Q^J^ The energies of the initial and final electron, 
of the tagged photon and of the initial proton (e, 62, and Ep) are defined in the labora- 
tory reference frame (i.e., the rest frame of HERA detectors). The cross section (j^) agrees 
with [^, ^. Also note that we explicitly included the correction from the vacuum polarization 
operator 11 (—Q^) in the virtual photon propagator. The aim of our work is to calculate the 
higher order QED radiative corrections for this process in the leading and next-to-leading 
logarithmic approximation. 

In this paper we restrict ourselves to the model independent QED radiative corrections 
related to the lepton line, which form a complete, gauge invariant subset for the neutral 
current scattering process. The remaining source of QED radiative corrections at the same 
order, such as virtual corrections with double photon exchange and bremsstrahlung off the 
partons are more involved and model dependent, they will be considered elsewhere. Our 
approach to the calculation of the QED corrections is based on the utilization of all essential 
Feynman diagrams that describe the observed cross-section in the framework of the used 
approximation. The same approach was used recently for the calculation of the QED correc- 



tions for the small angle Bhabha scattering cross-section at LEPl [IC]. This work extends 



our previous calculation at the leading logarithmic level |11| and presents the details of the 



brief outline published in |12 



The paper is organized as follows. Section |^ is devoted to the corrections related with 
emission of virtual and soft real photons in the hard coUinear photon emission DIS process. 
In sect. we consider the radiative corrections due to emission of two hard photons in the 
collinear kinematics (where we distinguish between the cases when both photons are emitted 
close to the initial electron direction and the case when one of the photons is emitted along 
the initial and the other one along the scattered electron direction) and the semi-coUinear 
kinematics, where the additional hard photon is emitted at a large angle. Section ^ collects 
the results obtained and discusses two experimental cases: an exclusive set-up, that assumes 
that a bare electron can be measured, and a calorimetric one. We show that, in the latter 
case for a coarse detector that performs a calorimetric measurement, we can reproduce the 
result of our previous paper , where the leading logarithmic approximation was used and 
where the emission along the final electron was not taken into account. In conclusion we 
give some numerical estimates. The appendices are devoted to details of the calculation. 



^The corresponding Born cross section including contributions from the Z can be found in ref. 
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2 Virtual and soft corrections 

In order to calculate the contributions from the virtual and soft photon emission corrections, 
we start from the expression for the Compton scattering tensor with a heavy photon ||T3[| , 

K,, = {8na)-' ^ Mf-^'^iM^^'-^'^y , (8) 

spins 

where is the matrix element of the process of Compton scattering 

7*(-g) +e(pi) ^7(A;) +e(p2), (9) 

and the index describes the polarization state of the virtual photon. This tensor is conve- 
niently decomposed as follows: 

K,. = 1{P,u + P:,), (10) 

Zn Zn 

Qfiu = g^lu Pul = Pifi - Qf^—^ , P2i, = P2i, - qt,—^ , pi = q + p2 + k . 



The expressions for the quantities Bij corresponding to the Born approximation are: 



2 

I - / 1 1 \ A ^2 Q^2 



(s + uY + {t + uY - 2m'q' ^ + , Bu = 



1 1\ ^ V 8m' 



2 tV St S2 



B22 = — s = 2p2-k, t = -2pi-k, u = {p2-piY, 

g2 = s + t + u, p\=pl = m^, k'^ = 0. (11) 

The one- loop QED corrections are contained in the quantities Tij, whose explicit expressions 
are given in |]13[. Here we have to integrate them over the solid angle of the emitted photon 
corresponding to the shape of the photon detector. We need to keep only the terms singular 
in the limit — >■ 0, since after integration the constant terms contribute only proportional 
to 6q ~ 10~^ and can be safely neglected. Another simplification comes from the fact that 
we need only the symmetric (and real) part of the tensor K. This way, by using typical 
integrals 

dQk 1 I/O f dQk m? 1 ^^^-j 



27r t 2e2(l-^)' J 27r 2e^{l-zf' 

and using the expressions given in Appendix A we obtain the following expression for the 
Compton tensor integrated over the angular part of the photon phase space: 



l + _,)p(,,i„)__T 



(13) 



^We have already dropped those terms that vanish in the high-energy Umit when one integrates over any 
finite region of photon phase space. 
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p = 41nA(LQ-l)-Lj + 3LQ + 31n^+^-^, 

777- O ^ 

1 -I- 4z 2 — fl — z")^ 

T = -^{Alnz + B)-- LqIuz-—^^ ^Lq + O(const) 

1 — z 1 — z 21 — 2; 



A = -Ll + 2LoLq - 2Lo\n{l - z), B = (\n^ z - 2U2{1 - z)j Lo 



U,{x) = - I ^ln(l-y). 



1^ 2 



The quantity A, which enters into the expression for p, is a fictitious photon mass. 

In the construction of the total expression for the tensor K^^ we replaced = qu = 
P2fj,,u = zpifj,^u, bearing in mind the gauge invariance of hadronic tensor [|14| , 



Ayr / ~ ~ \ 

= —{W2{xh,Ql)P,P.-M'W,{xH,Ql)g,,), 



H,, = ^(W2{x,,Ql)P,P,-M'W^{xh,Qlrg,u), = , (14) 



2P-qh 



J^u qhf 9 
1h 



Here we imply qh = q, Qh = —q^- 

Consider now the process with emission of a soft photon in addition to the emission 
of the hard one, which hits the PD. We imply the condidtion that the energy of the soft 
photon should be less than some small quantity 6e (in the center-of-mass system). In 
straightforward calculations, starting from Feynman diagrams, some care is to be paid in 
the evaluation of integrals over the phase volume of the soft photon, as some contributions 
are crucially dependent on the correlation between our two small parameters A = 6e/e and 
^0- In our particular case 6'o <C A -C 1, the result coincides with the one obtained using 
the approximation of classical currents for soft photons. The total effect for the sum of 
contributions of virtual and soft photon emission consists in the replacement of the quantity 
p by p in eq. (|T3|) (see eq. (45) in |[r3|): 



p ^ p = 2(Lq - 1) In^ + 3Lq + 31n^ - In^F - ^ - ^ + 2Li2(^) . (15) 

The final expression for the virtual and soft photon emission corrected tagged photon 
cross-section has the form 

^-f^=(^)'lP(z,U)^-T]t. (16) 

yaxayaz kit: J 



3 Double hard bremsstrahlung 

Consider now the emission of an extra photon with the energy more than 6e. For the 
calculation of the contributions from real hard bremsstrahlung, which in our case correspond 
to double photon emission with at least one photon seen in the forward detector, we specify 
three specific kinematical domains: i) both hard photons strike the forward photon detector, 
i.e., both are emitted within a narrow cone around the electron beam {9 < Oq); ii) one 
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hard photon is tagged by the PD, while the other is coUinear to the outgoing electron 
(6*2 = ^2^2 — ^0)5 finally iii) the second photon is emitted at large angles (i.e., outside 
the defined narrow cones) with respect to both incoming and outgoing electron momenta. 
We denominate the third kinematical domain as a semi-coUinear one. The contributions 
of the regions i) and ii) contain leading terms (quadratic in the large logarithms Lq, Lq), 
whereas region in) contains formally non- leading terms of order Lq ln(l/^Q), which, however, 
give a contribution numerically larger than the leading ones since e6o/m -C I/Oq. 

The calculation beyond the leading logarithmic approximation may be performed using 
the results of a paper of one of us |jl6|. The contribution from the kinematical region i) 



(both hard photons being tagged), has the form (see eq. (11 6) from 



3^77 



y dx dy dz 



a 



+ 



In^z 



- - 2 In A 
2 



4ii±il!ln^-^ 



))+6(l-z) + 



A 



+ C (const) 



(17) 



Here we use the notation Pq^ (z) for the 0-part of the second order term of the expansion 
of the electron non-singlet structure function 

5(1 -z) + ^P^'\z)L + I ( ^yp(^)(z)L' + ..., 



D{z,L) 
P^'\z) 



2tx ^ ' 2 V27r 
P^\z)Q{l~z-/\) + Pf'^5{l-z) 
l + z^ Ml) 3 



A ^ 0, 



1-z' 



1-z 



p> 



+ 21nA 



3\ 1 

21n(l - z) - hiz + - \ + -{I + z)\nz - I + z 



(18) 



The parameter A serves as the infrared regularization parameter. 

The contribution of the kinematical region ii) to the observed cross-section depends on 
the event selection; in other words, on the method of measurement of the scattered particles. 

In the case of exclusive event selection, when only the scattered electron is detected, 
while the photon that is emitted almost collinearly (i.e., within a small cone with opening 
angle 29'q around the momentum of the outgoing electron) goes unnoticed or is not taken 



into account in the determination of the kinematical variables, we have (see 118 from |]T6| ) 



z d'^a. 



3^77 



y dx dy dz 



a 

4^ 



P{z,Lo) 



dy2 



A/y 



I + 2/2L 



l + (l + y2) 



1/2 



-(Z-1) 



y2 



j:{xb,yb,Ql 



(19) 



where 



\n{ee'Jmf + 2\nY, 1/2 = y , ¥ = 82/6, 



Xb 



xyz{l + y2) 
-(l-l/)(l + l/2) 



Vb 



Z-(l-y)(l + l/2) 



2z-Y{l + c) 
Y{l + c) ' 

Ql = Q'zil + y2). 



(20) 
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More realistic (from the experimental point of view) is the calorimetric event selection, when 
only the sum of the energies of the outgoing electron and photon can be measured if the 
photon momentum lies inside the small cone with opening angle 26q along the direction of 
the final electron. In this case we find 



3^77 



dv: 



i'Lcal 



y dx dy dz 



a 

47r^ 



P(z,Ln 



■P(z,Ln 



dy2 



A/Y 



(1+^/2)3 



1 + (1 + ^/2) 



(L- 1) 21n 



Y 
A 



1/2 

3 
2 



iL-l) + y2 



E. 



(21) 



In the last equation we used the relation 

= 



(1 + 1/2)^ 



(22) 



which is valid for the calorimetric set-up. 

Consider at last the semi-coUinear region in). The relevant contribution may be calcu- 
lated using the quasireal electron method |T5|: 



z dV2 



3^77 

Hi 



y dx dy dz 



—r[z,Lo)— / — — — I 

27r IT J UJ2 Q% 

Bp,{zp,,p2,h)HPy{87r). 



(23) 



The quantity Bp^{zpi,p2, ki) is obtained from equation (|TT|), where it necessary to set m = 0. 
After some algebraic transformations we obtain 



P 



1 



F2{xsc, Ql) ^XscG - Vlx^.iz' + (1 - y f)V + (1 - y){zQ^ - s) 



z{zQ^~t) I -GFi(x,„QL) 

Q 



G = z'Q^ - 2st + Ql 



(24) 



)2 

' sc 



T/^ , 1 ^ oB~r ' s = 2p2-k2, t = -2zpi ■k2, Ql^ = zQ'^ - s - 1 . 
V[z + y - 1) - 2P ■ k2 



The angular integration in eq. (^) is to be performed over the whole phase space, excepting 
the small cones along directions of motion of the initial and scattered electrons that corre- 
spond to the kinematic regions i) or ii). The result (for the details see Appendix B) has the 
form 



z d-^cr; 



3^77 



y dx dy dz 



a \ 



2ii 



P(z,Ln 



dX2 + {Z - X2)\ 2(1 -c) 

7 ^ In ^2 + 

X2 Z[Z - X2> t^o 



+ 



dX2l + (l+t/2)% 2(1 -C) ^ 

In — -77^ — L, 



X2 



1 + 2/2 



t) 5 



(25) 
(26) 
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The logarithmic dependencies on the infrared regulator A and on the angles ^o, fully 
contained in the first two terms on the r.h.s., whereas the quantity Z represents an integral 
over the whole photon phase space of a well-behaved function, and it is free from coUinear 
and infrared singularities. Its explicit expression is given in Appendix B. 
The upper limits of the X2-integration in (^) read 

(27) 

and the arguments of are 

xt = ■ -, yt = , Qt = Q{z-X2). (28) 

Z - X2 + y - I Z - X2 

An explicit expression for Xm, which is relevant for the calculation of Z, is given in Ap- 
pendix B. 

The formulae given above (see eqs. (|^, (p^Bj) , (|l^, (|19|) or (^), and (^SD ) provide the 



complete answer for the leading and subleading contributions up to the second order of 
perturbation theory. The total sum of virtual, soft, and hard additional photons emission 
corrections to the radiative DIS cross section does not depend on the auxiliary parameter 
A = fe/e, as it should be. 



4 Results for different experimental situations 

The sum of the contributions of the leading and next-to- leading corrections at order a^, 
which are given exphcitly in expressions (0), (|T^, (|T^) or (PT|), and (P^, may be written 
in the form 



z dV 



(29) 



y dx dy dz \2'k, 

The first term Ej is independent of the experimental selection of the scattered electron and 
has the form 

"1 - 16^-^2 



-LlPi\z) + P{z,Lo) 



-Fln^F -2Li2(-2) + 2Li2 



+P{z, Lo) t In 



2(1 + ^2 

l + c\ 2(1 + ^)2 



+ 3-21nF + 



Az 



l + z' 



In^; -|- 



1 + ^2 



r "^0 



du 



u 



du 



[l + {l-u) 







u 



no 



+ P(^,Lo)Z, 



u 



ln(l - z) + 

' St 
,{\-u)t 

X2 



Z 



Uo 



2(1 + ^2- 

A- 



x\ 



\n'z 



(30) 



where Z is given in Appendix B and the remaining notations are as above (see (p!7|) , (0), 
and (^). 

The second term in (p9D, denoted Sj, however, does explicitly depend on the event 
selection. It corresponds to the emission of a hard photon by the scattered electron. In the 
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exclusive set-up, when only the scattered bare electron is measured, while the photon that 
is emitted close to the final electron's direction is ignored, this contribution reads 



excl 
f 



ilQ + lnY-l] 



xllY 

Piz,Lo) J dy2 



xe (^2 - ^) + (Lq + \nY- 1)5(2/2) (2 In ^ + ^ 



1 + (1 + 1/2) 

y2 



y2 



1 



1+1/2 



(31) 



In this case the parameter 6'q, that separated the kinematic regions ii) and in), only plays 
the role of an auxiliary one; it has already cancelled in the above expression for the cross 
section. 

As we will see below, this situation is quite different for the experimentally more realistic, 
calorimetric set-up, when the detector cannot distinguish between events with a bare electron 
and events when the electron is accompanied by a hard photon emitted within a small cone 
with opening angle 26q around the direction of the scattered electron. For this case we obtain 



2 9'i I 2/2 1 + 1/2 



X s,e(^/^^-2/2)- 



,l+l/2) 



(32) 



For the calorimetric event selection the parameter 6'q is a physical one and the final re- 
sult therefore does depend on it. However, the mass singularity that is connected with 
the emission of the photon off the scattered electron is cancelled in accordance with the 



Kinoshita-Lee-Nauenberg theorem 17 



Note that the case of a coarse detector for the scattered electron, i.e., 6*0 ~ 



at the level of leading logarithms with the result of paper |11 



(9(1), agrees 
that was obtained in the 
approximation of absence of emission along the scattered electron. Our result disagrees with 
the result of Bardin et al. [0] on the radiative corrections, as they neglected the interference 



of the emission of two photons; see |Tl| for a detailed discussion. 

Finally we will give some numerical values obtained for the radiative corrections at lead- 
ing and next-to-leading order for the experimentally relevant case of the calorimetric event 
selection with a realistic resolution. As input we used 



E, = 27.5 GeV . 



Ep = 820 GeV , Bo = 0.5 mrad , 



(33) 



whereas for the resolution of the detector we assumed 6'^ = 50 mrad. As structure function 



we chose the ALLM97 parameterization with R = 0. No cuts were applied to the photon 
phase space. 

Figure |l| compares the radiative correction 



dV 



JRC 



dVeo 



- 1 



(34) 



at leading and next-to-leading order for the measurement in terms of the standard Bjorken 
variables x and y for x = 0.1 and x = 10~'^ for a tagged energy of EpD = 5 GeV. The typical 
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size of the next-to-leading order contributions with respect to the leading-logarithmic result 
]TT| amount to up to the order of 5%, depending on the kinematic region. The apparent 
cutoff in the curves at low y appears at 



1/Bj,min = (35) 

1 — xz 

due to the condition that a; < 1, see eq. (|^). 

For the case of the shifted Bjorken variables (|^) we present the corresponding result 
in fig. H. Here no lower cut in y appears, like in the case of standard Bjorken variables. 
However, the curves for small x = 10~^ are not continued below the value of y where the 
cones with opening angles 26*0 (PD) and 26*0 (calorimeter resolution) would start to overlap, 
as y — ^ corresponds to forward scattering. The difference of the radiative corrections 
between leading and next-to-leading order in this case also amounts to the order of 5%. 

Lastly, in order to exhibit the 2;-dependence of the next-to-leading order contributions, 
i.e., the dependence on the energy in the photon tagger, we plot the corrections for E'pd = 
20 GeV in fig. |^ for x = 0.1 and x = 10~^. One sees in particular the increasing relevance of 
the next-to-leading order terms for large Epi) and in the experimentally interesting range of 
small X. 

We note in conclusion that the set of Feynman diagrams considered here is gauge invariant 
and model independent but not complete. We have neglected the contributions with two 
virtual photons exchanged between electron and the target that appear at the same order 
of perturbation theory, as well as the interference with the contributions when the second 
photon is emitted by the hadronic side. (These contributions have not yet been considered 
in any hterature about DIS known to the authors). However, the description of this part is 
definitely model dependent and will be discussed elsewhere. 
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Appendix A 

List of angular integrals for virtual corrections 

In this section we collect the results of the angular integration of the definite structures of 
the Compton tensor in [0. Using integrals similarly to (0) and retaining only terms that 
contain at least one large logarithm Lq or Lq, we obtain 



2£2 

Of 



271 ' 



1 + 



-{Lq 



1) + 1 



l + z^ 
(1^ 



:[A\yiz + B] 



Az 



(1 



Lq In 2; 



2{l-zf 



Lq . 
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2e' 



27r 



Az 



11 



2^(1 + (1-^)2)^^, ,2(3-2) 
-pLo ^-—^—^{A\nz + B)- ^ + 



l-z)2""" (1-2^4 
2Lo fz{8z-3) 



[l-zY 



2e' 



'2^ 



22 



1 - 1 - -2 

Lq — 



^ (1-2)2 

2^(1 + 2(1 -z^2^ 



\nz + 2z + z"^ 
16 

+ 



'1-2) 



In zLo - 



(36) 



[1-zr 

2Ln /I +42(^2 + 2-1 



(Alnz + 5) - ^ + 



26^ 



z{l-z)^ 
2^2 



1-z 



z(l -z)= 
In 2; + - + 42; - 1 



21 



{1-zy 



3^ — 1 

{Alnz + B) + ^ -A + 



2e' 



~2^ 



(1-2)3 

2z(2-z 



{\-zY 

-1 +4:Z - 4:Z^ - AZ^ 



12 



\-z 

{AX^z^B) + - 



In ^ - 2;z^ - 2z + 1 



3 - 2 , 2Lo 
■A + 



3-8z 



;i-;z)4 ^ ~ ' ' ■ (1-^)3 V 1-z 

where p, A and i? are given by eq. (plsl). It is remarkable to see that the relation 



ln;z - 1 - 2z 



2ti 



AzTg + Q/(Tn + z'Tss + zT^2 + 2T21; 







(37) 



is fulfilled, leading to the factorization of the virtual corrections in eq. (|T^). 

Appendix B 

Angular integration for bremsstrahlung corrections 

To perform the angular integration in ( pBj) we first represent the integrand in the form 



Qi 



ut 



(3^ 



l''2 



where ti^2 = (1 — ci,2)/2, ci,2 = cos 6*1,2, and Oi^2 are the angles between non-tagged photon 
momentum k2 and the momenta of the initial and the scattered electrons. Note that F{ti, ^2) 
behaves regularly for ti — > or t2 0. This can be easily seen by considering the limiting 
cases for the quantity P . 

For the case t — >■ 0, which corresponds to the the second photon being emitted close to 
the direction of the incoming electron, one obtains from eq. (El) 



P\t^O = ^(Z' + (Z-X2f) 



X2t 



X 



xtF2{xt,Q] 



Ql x^y'^{z - X2) 



Fi{xt,Q 



(39) 
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while for the case s — 0, corresponding to the second photon being almost coUinear to the 
final electron, 



|s-»0 



y2s 



(1 + (1+^/2)') X 



XbF2{xb,Qt){--o - 



Ql x'^y'^z{l + y2) 



Fiixb,Q 



see eqs. fpoD and ( P8|) for the notation. The r.h.s. of eq. (|38D is easily seen to be 



Q 



1 a z'^ + {z - X2f 
tit2 IGvrxl z{z — X2) 

J g 1 + (1 + 1/2)' 

tit2 IQtTxI 1+1/2 



j:{xt,yt,Qt) 



j:{xb,yb,Ql) , 



where a = (1 — cos6')/2. 

For the phase space of the photon we use the following representation: 



J^=e'Jx2 dX2 dfi2 = 4.^ / X2 dX2 J ^©l^) 



dti dt2 



D = it2- y-){y+ - t2), y± =h{l-2a) + a± 2^a{l - a)U{l - h 
The region of integration is determined by the conditions 

CTi < ti < 1, CT2 < t2 < 1, D> 0, ai = a2 = 



Using the substitution 



t2 t2{ti,u) 



y+ + u^y^ 



and the identity 

1 1 



(Tl (72 



tlt2VD 



TT 

a 

+2 



F(a,0) In— + F(0,a) In — 

0"2 CTi 



+ 



d-u 



lim 

1 + u-^ v-^o 



dti 



(40) 

(41) 
(42) 



(43) 



(44) 
(45) 



ti\ti — a 



{Fiti,t2) - F{a,0)) + 



/Af 
^(F(a,0)-F(0,a)) 



which is valid for ai, o"2 ^ 0,, we obtain for Z from eq. ( pSD the following expression: 



(46) 



4(1 - c) 



/ ^ f / / ^(*(t. t2(t„ .)) - $(a, 0)) + 



+ /^/^($(a,0)-$(0,a)) 

J It a J X9 



tia J X2 

r? 



(47) 
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where 



$(tl,t2) 



Qtc 



ci=l-2ii, C2=l-2t2(ti,M), c=l-2a 



(48) 



The upper hmit of the X2-integration, Xm, may be deduced from |0. It has the form 



This finally leads to eq. 

It is important to note that when calculating Z one encounters neither coUinear nor 
infrared singularities. 
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